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The oxidation kinetics of CVD SiC were monitored by ther- 
mogravimetric analysis (TGA) in a 50% H 2 O/50% O z gas 
mixture flowing at 4.4 cm/s for temperatures between 1200° 
and 1400°C. Paralinear weight change kinetics were 
observed as the water vapor oxidized the SiC and simulta- 
neously volatilized the silica scale. The long-term degrada- 
tion rate of SiC is determined by the volatility of the silica 
scale. Rapid SiC surface recession rates were estimated 
from these data for actual aircraft engine combustor 
conditions. 

I. Introduction 


The rate of reactions (1) and (2) can be described by a para- 
bolic rate constant, k p . The rate of reactions (3)-(5) are given 
by linear volatilization constants, k { . When oxidation and vola- 
tilization occur simultaneously, the oxidation rate is affected by 
the volatilization rate, and the overall kinetics are described by 
the more complex paralinear kinetics. 6 

Paralinear oxidation kinetics can be monitored in terms of 
either changes in oxide thickness or changes in overall sample 
weight. Paralinear oxide scale thickness kinetics have been 
mathematically described by Tedmon for Cr 2 0 3 -forming Fe-Cr 
alloys. 7 This expression is directly applicable to SiC and is 
as follows: 


S ilicon CARBiDE-based composites are currently being devel- 
oped for use in the complex environments of heat engines. 
Possible degradation mechanisms for SiC in combustion envi- 
ronments have been reviewed by Jacobson. 1 While the behavior 
of SiC in oxygen has been well studied, less information on the 
durability of SiC in more complex combustion environments is 
available. Some burner rig studies of SiC degradation have been 
conducted; 2 " 4 however, specific degradation mechanisms are 
difficult to identify in this type of study since the environmental 
variables are not well defined and degradation rates are difficult 
to quantify. Water vapor is a major component of all combus- 
tion environments. Jacobson 1 has calculated that the products 
of combustion contain about 10% water vapor independent of 
the fuel-to-air ratio. In this study, we have attempted to isolate 
the basic mechanisms of SiC degradation in water vapor. 

In an environment containing water vapor and oxygen, ther- 
modynamic calculations indicate the primary reactions which 
occur for the oxidation of SiC are 


SiC + §0 2 (g) = Si0 2 + CO (g) 

0) 

SiC + 3H 2 0(g) = Si0 2 + 3H 2 (g) + CO(g) 

(2) 

We report in this paper the first observation that the silica scale 
formed on SiC in an oxygen/water vapor gas mixture may 
simultaneously volatilize by forming a silicon hydroxide or 
silicon oxyhydroxide species under atmospheric pressures. 
Some possible volatilization reactions are 5 

SiO, + H 2 0(g) = SiO(OH) 2 (g) 

(3) 

Si0 2 + 2H 2 0 (s) = Si(OH)„(g) 

(4) 

2Si0 2 + 3H 2 0(£) = Si 2 0(0H) 6 (g) 

(5) 


Silica volatility in water vapor will be discussed in more detail 
below. 


Nathan S. Jacobson — contributing editor 



where t is the oxidation time, jc is the oxide thickness, k' p is 
the parabolic rate constant for oxidation in units of thickness 2 
per unit time, and k[ is the linear rate constant for scale volatil- 
ization in units of thickness/time. At long times, the scale 
thickness becomes a constant as the oxidation rate equals the 
volatilization rate. This limiting scale thickness, x f , is equal to 
k’/2k[. The oxide thickness as a function of time as predicted 
by Eq. (6) is shown in Fig. 1 for typical values of k’ v and k[. 

Similarly, paralinear kinetics for oxidation of SiC can be 
expressed in terms of measured weight changes 6 using the 
following two expressions: 


cc%\ 2*, Aw, _ 2 kAwX 

' - wr^r ~ "r ^ 

Aw 2 = — fU,/ 


(7) 

(8) 
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Fig. 1 . Paralinear oxide growth kinetics typical for SiC in 50% water 
vapor flowing at 4.4 cm/s and temperatures between 1 200° and 1 400°C. 
The curve is a graphical representation of Eq. (6). After about 3000 h, 
a limiting oxide thickness equal to k f p /2kl is achieved. 
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where A w x is the weight gain due to growth of the scale, A w 2 is 
the weight loss of the scale due to volatilization, k p is the 
parabolic rate constant for oxidation in units of mass 2 /area 2 
time, k x is the linear rate constant for SiO z volatilization in 
units of mass/area time, a = MW Si02 /(MW Oi - MW C ), and 
P = MW SjC /MW Si02 where MW is the molecular weight of the 
species shown in subscript. The total observed weight change 
is the sum of Aw^ and Aw 2 . At long times when the limiting 
oxide scale thickness is achieved, the weight change is simply 
expressed by Eq. (8). The weight changes for SiC in 50% 
H 2 O/50% 0 2 predicted by Eqs. (7) and (8), as well as their sum, 
are shown in Fig. 2 for typical values of k p and k x . These rate 
constants, though expressed in terms of weight change rather 
than oxide thickness, are equivalent to those used in Fig. 1. 
Since the abscissa is the same in Figs. 1 and 2, direct compari- 
sons of oxide thickness and weight change at any given time 
can be made. Note that the weight loss is nearly linear at times 
well before scale growth and vaporization rates balance one 
another and the limiting oxide thickness is achieved. 

Some information on both the parabolic and linear rate con- 
stants is found in the literature and will now be summarized. 
Hie parabolic oxidation rate of silicon carbide in oxygen by 
reaction (1) has been extensively studied. 8,9 The parabolic oxi- 
dation of silicon carbide in water vapor has also been studied 
but is somewhat less understood. A summary of these studies is 
given in Ref. 10. While these studies in general show the rate of 
oxidation is enhanced in water vapor, a systematic study of the 
dependence of the oxidation rate on water vapor partial pressure 
has not yet been completed. However, it has been shown that in 
H 2 0/0 2 gas mixtures, water vapor is the primary oxidant for 
Si. 11 It has also been shown that the parabolic oxidation rate is 
proportional to the water vapor partial pressure for Si, 11 and 
preliminary results indicate the same holds true for SiC. 12 

The oxidation/volatilization of SiC has been reported by Kim 
and Readey 13 and Narushima et ai ]A in cases where H 2 /H 2 0 
or C0/C0 2 gas mixtures are used. In these cases the oxide is 
reduced by hydrogen or CO to form the volatile SiO species. In 
the work by Kim and Readey, the minimum hydrogen partial 
pressure was 0.9 atm. It must be emphasized that in our work 
the maximum hydrogen partial pressure from the decomposi- 
tion of water was 10~ 5 atm. Thermodynamic calculations show 
that SiO is not present in large enough quantities to account for 
the observed silica volatility in the oxidizing conditions studied 
here; therefore, some other volatile species must account for the 
observed paralinear kinetics. 

The volatility of free-standing silica in high pressure (on the 
order of 100 atm) water vapor has been known for some time 
and is summarized by Glemser and Wendlandt. 15 It has been 
proposed that volatilization of silica under these conditions 



Time (h) 

Fig. 2. Paralinear weight change kinetics typical for SiC in 50% 
water vapor flowing at 4.4 cm/s and temperatures between 1 200° and 
1 400°C. The dotted curves are from Eqs. (7) and (8), while the solid 
curve represents the total weight change expected. 
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occurs by reactions (4) and (5). Other work by Cheng and 
Cutler 16 and Hashimoto 17 has shown that silica is also volatile in 
atmospheric pressures of water vapor. The careful transpiration 
study by Hashimoto 17 shows a square dependence of the volatil- 
ity rate on water vapor pressure, indicating that reaction (4) is 
dominant. The likely volatile species under these conditions is 
thus Si(OH) 4 . Preliminary mass spectrometric results at atmo- 
spheric pressures confirm the presence of this vapor species. 18 
Finally at very low water vapor pressures, 10“ 5 atm, a Knudsen 
cell mass spectrometric study by Hildenbrand and Lau 19 has 
identified SiO(OH) 2 as the volatile species, indicating that 
reaction (3) may be important under these conditions. 

The idea of simultaneous volatilization of silica formed dur- 
ing oxidation of silica-forming materials (Si, SiC, Si 3 N 4 ) in 
water vapor/oxygen gas mixtures follows naturally from the 
above volatility studies of free-standing silica, yet there are 
very few observations of this in the literature. The volatility of 
silica formed as an oxidation product on silicon above 25 -atm 
water vapor pressure has been noted in the literature. 20 Some 
evidence of the erosion of silica thermally grown on silicon in 
1 atm wet oxygen (83% H 2 0) and steam was observed. 21 

The objectives of this paper are, first, to demonstrate conclu- 
sively that paralinear oxidation of SiC does occur in water 
vapor. Once this has been demonstrated, the second major 
objective of this work is to accurately determine the silica 
linear volatilization rate constant in water vapor. Finally, the 
volatilization rate will be used to predict the long-term durabil- 
ity of SiC in combustion environments. 

II. Experimental Procedure 

( 1) Material 

The material used in this study was polycrystalline mono- 
lithic CVD (3-SiC (Raytheon/Bomass Machine Specialties, 
Somerville, MA). The major impurities as determined by semi- 
quantitative dc arc spectrographic emission were Be and Cu, 
both in quantities of less than 100 ppm. The sample coupons 
had dimensions of approximately 2.5 X 0.7 X 0.3 cm and had 
a 600 grit surface finish. Some experiments were also con- 
ducted with fused quartz coupons (Quartz Scientific, Fairport, 
OH). Manufacturer’s specifications state that this material has a 
nominal purity of 99.995%. These sample coupons had dimen- 
sions of approximately 2.5 X 1.2 X 0.1 cm. All samples were 
ultrasonically cleaned in detergent, deionized water, acetone, 
and alcohol prior to furnace exposures. 

(2) Thermogravimetric Experiments (TGA ) 

The durability of SiC in water vapor was studied between 
1200° and 1400°C. Initial studies of SiC had been conducted in 
a 10% H 2 0/balance 0 2 gas mixture 10 which is close to the water 
vapor level found in the combustion products at 1 atm. Little 
effect of H 2 0 vapor at this level was found on the oxidation of 
SiC. It was thought that a 50% H 2 O/50% 0 2 gas mixture at 
1 atm total pressure would more realistically simulate combus- 
tion conditions which are typically at higher pressures (e.g., 
10% H 2 0 at 5 atm total pressure). At the low flow rates typically 
used in TG studies, 0.4 cm/s, the water vapor saturation level in 
the gas stream was found to be less than 50% even with the 
water bath temperature at the desired dewpoint of 81.7°C. Sur- 
prisingly, at a constant bath temperature as the flow rate was 
increased, the saturation level also increased. A total flow rate 
of 4.4 cm/s was used to obtain 50% saturation of the oxygen 
stream. The highest water vapor partial pressure that was repro- 
ducibly obtained using the water saturator already described 10 
was 0.5 atm. Fortuitously enough, use of a higher flow rate led 
to the observation of silica volatility. 

Three types of thermogravimetric experiments were con- 
ducted to determine the parabolic oxidation rate ( k p ) and linear 
volatilization rate (k x ) constants, as explained below. In all 
cases, samples were weighed before and after each TGA experi- 
ment. First, a cleaned as-received SiC coupon was suspended 
from a microbalance (Cahn 1000, Cahn Instruments, Cerritos, 
CA) using a sapphire hanger. The samples were suspended 
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within a 2.2 cm ID fused quartz reaction tube (Quartz Scien- 
tific). Fused quartz reaction tubes were used, since previous 
work had shown that the use of high-purity alumina tubes with 
water vapor substantially increased A p for SiC, 10 due to residual 
volatile Na impurities. The samples were heated using a MoSi 2 
furnace. The hot zone of this furnace was about 4 in. 
(~10.2 cm) long. The sample temperature was calibrated to 
the furnace setpoint in separate experiments, using a type R 
thermocouple. The thermocouple was removed during experi- 
ments with sample coupons to avoid contamination from the 
alumina thermocouple sheath. The temperature uncertainty is 
±5°C. A gas mixture of 50% O 2 /50% H 2 0 flowed through the 
reaction tube at a room temperature rate of 1 L/min (4.4 cm/s). 
A counter flow of 0.4 L/min dry argon was circulated through 
the balance chamber to prevent water condensation in the 
microbalance. Gas flow was controlled using a Tylan mass flow 
system. Water content of the gas stream was monitored by 
measuring the condensate in the gas exhaust daily and was very 
constant. Weight and time were recorded every 6 min for 100 h, 
using a data acquisition system. Equations (7) and (8) were 
fit to the resulting specific weight-vs-time curve, using the 
procedure described below. The rate constants A p and A, were 
obtained from the best-fit curve. 

A second experiment was conducted to directly obtain the 
linear volatilization rate of silica, A,, under the same conditions. 
A fused quartz coupon was exposed to the water vapor/oxygen 
gas mixture, using the same TGA procedure as described in the 
previous paragraph. A linear weight loss uncomplicated by the 
oxidation reaction for SiC was obtained. The rate constant, A,, 
was determined from the final weight minus the initial weight 
and was in good agreement with the rate determined from a 
linear regression of the weight loss data. 

A third type of thermogravimetric experiment was also con- 
ducted to determine the linear volatilization rate constant. Here, 
the SiC was preoxidized to grow a thick oxide scale (6 to 8 pm). 
Ideally, a preoxidation exposure in dry oxygen in the TGA 
would be the simplest procedure. However, the time required to 
grow a thick oxide in this manner is inordinately long. Since 
water vapor enhances the oxidation rate of SiC, it was found 
that a thick oxide scale could be grown more rapidly in a 
horizontal tube furnace in 50% H 2 O/50% 0 2 in which condi- 
tions were such that the gas stream was saturated with the 
silicon hydroxide species, thus suppressing the silica volatility. 
This preoxidation procedure was conducted simultaneously for 
four coupons and is described in detail elsewhere. 12 After pre- 
oxidation, an approximate oxide thickness for each sample was 
determined from the weight change. One sample from each 
preoxidation batch was cross-sectioned to measure the oxide 
thickness by SEM. No oxide spalling was observed. Each intact 
sample was inserted into the TGA and exposed to the 50% 
water vapor mixture. During this exposure to water vapor, the 
weight gain due to oxidation was negligible compared to the 
linear weight loss of the silica scale. The rate constant, A,, was 
determined from a linear regression of the weight loss data. 

(3) TGA Control Experiments 

Fused quartz coupons were exposed to 100% dry 0 2 in the 
TGA under identical conditions (1200°-1400°C, 4.4 cm/s total 
flow) to ensure that any volatility effects were due to the 
water vapor effects only. No weight change was observed in 
these experiments. 

There was concern that volatility of the fused quartz reaction 
tube would affect the oxidation kinetics of the SiC. Two specific 
issues were of concern. First, there was evidence that silica was 
depositing from the vapor phase on the sapphire hanger during 
the TGA exposure. Control experiments in 50% water vapor/ 
50% oxygen mixtures were therefore conducted. Sapphire 
hangers were suspended in the hot zone without samples and 
exposed to the water vapor/oxygen mixture at temperature. The 
TGA results show that linear weight gains were observed as 
volatiles from the fused quartz tube deposited on the portion of 
the hanger that was relatively cool. These deposits were identi- 
fied as silica, based on EDS results which showed silicon and 


oxygen. The linear weight gain rate was (1.37 ± 0.55) X 
10 2 mg/h* for six control experiments independent of tempera- 
ture between 1200° and 1400°C. The weight gains of the control 
runs were also in agreement with the difference between the 
TGA Awt and the Awt measured from the coupon before and 
after exposure, which was assumed to be linear and averaged 
(1.74 ± 0.33) X 10 -2 mg/h* for 19 experiments of all experi- 
ment types involving test coupons. Linear weight gains based 
on the weight difference technique were subtracted from the 
raw TGA data obtained for each coupon. 

The second concern was that the volatility of the fused quartz 
tube would affect the volatility of the test coupon. If Si-OH 
vapor species from the fused quartz tube saturated the gas 
stream before reaching the test coupon, volatility of the silica 
forming on the SiC coupon would be partially or completely 
suppressed. Therefore, the test coupon was placed near the 
upstream edge of the hot zone so that most of the volatile area 
of the tube would be downstream of the sample. The gas would 
have little time to be saturated with Si-OH species before 
reaching the test coupon. 

(4) Characterization Techniques for Oxidized Coupons 

Oxide thicknesses for each coupon were measured in cross 
section by scanning electron microscopy (SEM). Oxides were 
polished off one edge of each coupon using a 15 diamond 
wheel. Oxide thicknesses were measured approximately every 
millimeter along the edge of the coupon, giving about 25 total 
thickness measurements for each coupon. The accuracy of these 
thickness measurements was within 7% of the true thickness, 
based on calibration with an NBS SEM magnification standard. 
Some test coupons were also examined for surface morphology 
by SEM and phase composition by X-ray diffraction analysis 
(XRD). 

(5) Paralinear Weight Change Data A nalysis Methods 

The kinetic rate constants k p and A, for each paralinear data 
set were determined by nonlinear least-squares analysis and the 
paralinear kinetic model described in the Introduction. The 
objective function, x, to be minimized was 22 


x 2 (M.) 


^ A< bs - Aw ca1c (/'; A P ,A,) 
,-=i cr l 


(9) 


where Aw° bs is the experimentally observed specific weight 
change measured at time t t during the TGA run. The quantity 
Aw calc (f';A p ,A,) is the specific weight change calculated from 
the paralinear kinetic model for time t,. The variance, a,, was 
assigned a value of unity, since the true values associated with 
the experimental data were not known a priori. 

The quantity Aw ca,c (t,\ A p ,A,) is given by the relationship 

A vt ; j” alc = Aw,. + Aw 2[ (10) 


The quantities Aw, and Aw 2 are given by Eqs. (7) and (8), 
respectively. The differential form of Eq. (7) was used to com- 
pute Aw, by approaching the solution as an initial value prob- 
lem and then applying a 4th order Runge-Kutta integration 
technique 22 to solve for Aw, at time t { . Minimization of the 
objective function, x 2 (^ P ^\) was accomplished by simultane- 
ous adjustment of A p and A, using the Levenberg-Marquardt 
technique. 22 

The statistical uncertainties associated with the A p and A, 
estimates were determined by the “bootstrap” method, 22 23 a 
Monte Carlo-based simulation technique. In this method, the 
probability distributions associated with the estimated kinetic 
parameters may be reproduced from the original experimentally 
obtained data set by randomly resampling the original data set 
along with data point replacement. 23 This technique has been 
successfully used in many scientific disciplines ranging from 
astrophysics 24 to elementary particle physics, 25 as well as mate- 
rials science. 26 


‘Uncertainties reported are lo. 
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The generation of synthetic data sets in the bootstrap method 
consists of randomly selecting N data points from the original 
data set, R {0) = [p and then placing each randomly 
selected point into the N - membered set, R (l) . Once the random 
selection is complete, R {1) is now a subset of R m containing 
members of R {0) which are replicated with a probability of **\/e 
in R {[) . Therefore, R (]) will contain some members of R (0) which 
now will appear more than once while other members of R (0) 
will be absent in /? (1) . This process is repeated M times to 
produce M synthetic data sets. 

Each of the synthetic data sets were subjected to nonlinear 
least-squares analysis where M values of &p B) and k\ B) estimates 
were obtained. The mean and standard uncertainties were com- 
puted from the bootstrap estimates by the respective relations 


k 


(p.i> 


1 M 

_L V/. 

Mpr (pj) ‘ 


do 


and 


>-*■> 


M - 1 


M 


~ *(p.d) 


( 12 ) 


For typical production runs, M = 200. However, to determine if 
the estimated uncertainties determined by the bootstrap method 
were normally distributed, runs with M = 400 were performed. 

( 6) Calculations of Flux of Volatile Species 

The flux of the volatile species under different conditions can 
be calculated as follows. The flux of a volatile species through 
a gaseous boundary layer under laminar flow conditions over a 
flat plate, /, can be calculated using the following relationship: 27 



(13) 


where v is the linear gas velocity, L is the sample length parallel 
to the flow direction, is the gas viscosity, D is the interdiffu- 
sion coefficient for the diffusing species in the major gas spe- 
cies, p is the concentration of the diffusing gas species at the 
solid gas interface, and p' is the concentration of the major gas 
species. The gas concentrations are calculated using the ideal 
gas law. The interdiffusion coefficient was calculated using 
Chapman— Enskog theory. 28 Gas viscosities and Leonard-Jones 
parameters were obtained from the compilation of Svehla. 29 
Since these parameters are not available for Si(OH) 4 , values for 
SiF 4 were used as an approximation instead. 

Under conditions where gas velocities are high, turbulent 
flow occurs rather than laminar flow. The boundary layer is 
swept away and increased volatility rates are observed. In the 
extreme case where no boundary layer is formed, the maximum 
mass flux of the volatile species, /, due to free evaporation into 
a vacuum is given by the Langmuir equation: 


r „ / M 

y_/ V2i jRT (14) 

where P is the equilibrium vapor pressure of the volatile spe- 
cies, M is the molecular weight of the volatile species, R is the 
gas constant, and T is the absolute temperature. Equations (13) 
and (14) thus describe the minimum and maximum possible 
flux of volatile species. 


III. Results 

(1) Weight Gain and Scale Growth 

A summary of the weight gain (measured before and after 
TGA experiments) and oxide thickness results for the oxidation 
of as-received SiC coupons is shown in Table I. There is poor 
reproducibility for both the weight gain and oxide thickness at 
each temperature. Experimental variations of sample proper- 
ties, temperature, gas flow rate, and water content of the gas 
stream cannot explain this scatter. The high gas flow rates 


(4.4 cm/s) increased the amplitude of the noise in the TGA 
traces but did not affect the accuracy of the measured weight 
changes. Two possible contributions to this scatter are as fol- 
lows. Subtraction of the linear rate of silica deposition on the 
hanger from the overall weight change is a source of error in 
the calculated rate constants. The magnitude of this uncertainty 
has already been reported in the section on TGA control experi- 
ments. Another source of variability in the volatility rate among 
the sample coupons under otherwise identical conditions may 
be the differences in position in the hangdown tube. Positions 
off of true center of the small-diameter reaction tubes may 
result in interactions between the gas boundary layers of the 
sample and the tube wall. This would result in variations in the 
gas velocity and ultimately variations in the volatility rate from 
experiment to experiment. 

Besides the scatter in the weight gain and oxide thickness 
among the coupons, there was also a clear discrepancy between 
the weight gain and oxide thickness for any given coupon, as 
shown in Table I. For oxidation without volatilization, the 
weight gain and oxide thickness can be converted using the 
expressions 

Awt X 12.9 = x for cristobalite (15) 

Awt X 13.6 = x for amorphous silica (16) 

where Awt is in mg/cm 2 and x is in pm. These expressions 
already account for weight loss due to CO formation during 
oxidation. The uncertainty of thickness results based on SEM 
measurements will typically be ±10%. The variation in oxide 
thickness along the length of each coupon is greater than the 
uncertainty due to the technique. In this study, at 50% water 
vapor it was found that the ratio (measured Awt/predicted Awt 
from the measured oxide cross sections) varied between 0.09 
and 0.87, as shown in Table I. In all cases, the measured oxide 
thickness was greater than that predicted from the weight gains. 
In contrast, in previous work conducted in dry oxygen and 10% 
water vapor, 10 it was found that the ratio of the measured Awt 
to predicted Awt from the above expressions and the measured 
oxide cross sections varied from 0.96 to 1 .27. These data have 
not been reported before and are summarized in Table II. In 
10% H 2 0 and in dry oxygen, the agreement was good between 
measured weight gains and the measured oxide thicknesses. If 
anything, converting from scale thickness to expected weight 
changes slightly underpredicted the weight gain. 

The values of k p in 50% water vapor obtained from the fit of 
the paralinear data to the model are listed in Table III. These 
values are in fair agreement with values obtained by other 
methods 1012 that are available at this same water vapor partial 
pressure, as shown in Fig. 3. Values extrapolated from Ref. 10 
are obtained by multiplying by five the results which were 
obtained in 10% water vapor. This is done because preliminary 
results show that k p is approximately proportional to / > (H 2 0). 12 
The values from Ref. 12 were obtained in 50% water vapor 
under conditions where silica volatility did not occur. In this 
study, an activation energy for the oxidation of SiC in water 
vapor of 35 ± 147 kJ/mol (95% confidence interval) was 
obtained. Though the uncertainty is large, this result is not 
statistically significantly different from the value of 41 kJ/mol 
obtained in Ref. 10, the value of 62 kJ/mol obtained in Ref. 12, 
and the value of 68 kJ/mol obtained by Deal and Grove for the 
oxidation of silicon in water vapor. 1 1 

(2) Weight Loss and Scale Volatility 

The weight loss rates for both the fused quartz coupons and 
preoxidized SiC coupons were more reproducible than those 
obtained from the fit of the paralinear weight change of 
as-received SiC coupons. Linear volatilization rate constants 
were taken from the final minus initial weight for the fused 
quartz coupons and from the linear regression of the weight loss 
data for the preoxidized SiC coupons. The weight loss for the 
first 20 to 30 h was typically not linear in the experiments 
involving the preoxidized SiC coupons and was not included in 
the linear regression. 
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Table I. Comparison of the Total Weight Change, Awt, and the Oxide Thickness, x , for the 
Paralinear Oxidation of SiC in 50% H 2 O/50% 0 2 at 4.4 cm/s* 


Temp 

CO 

Time 

<h) 

Awt 

(mg/cm 2 ) 

,r f 

(txm) 

Awt 
from x 

Awt exp/Awt 
from x 

1200 

100 

0.113 

2.52 ± 0.45 

0.185 

0.61 


100 

0.226 

4.21 ± 0.36 

0.309 

0.73 


100 

0.214 

3.90 ± 0.57 

0.286 

0.75 


100 

0.192* 

Missing 




115 

0.018 

2.72 ± 0.55 

0.199 

0.09 


98 

0.059* 

2.98 ± 0.53 

0.219 

0.27 

1300 

100 

0.229* 

4.23 ± 0.65 

0.327 

0.70 


100 

0.152* 

3.68 ± 0.47 

0.284 

0.54 


100 

0.205* 

3.04 ± 0.54 

0.235 

0.87 

1400 

100 

0.243 

4.50 ± 0.67 

0.349 

0.70 


100 

0.271 

4.42 ± 0.54 

0.342 

0.79 


100 

0.050* 

3.31 ± 0.73 

0.256 

0.20 


100 

0.081* 

4.36 ± 1.02 

0.337 

0.24 


♦Conversion from x to Awt assumes density of amorphous silica at 1200°C (2.2 g/cm v ) and density of cristobalite at 1300 and 
1400°C (2.32 g/cm’). Uncertainties reported for oxide thicknesses are la for at least 20 measurements at 1 mm spacing along the 
SiC/oxide cross section, indicates samples which correspond to paralinear fits listed in Tables III and IV. 


Table II. Comparison of the Total Weight Change, Awt, and the Oxide Thickness, x, for the 
Parabolic Oxidation of SiC in 10% H 2 O/90% 0 2 and Dry Q 2 at 0.4 cm/s* f 


Temp 

(°Q 

Time 

(h) 

Awt 

(mg/cm 2 ) 

X* 

(pm) 

Awt 
from v 

Awt exp/Awt 
from x 

1200 

97 

0.164 

10% H-,0/90% 0 2 

1.90 ±0.23 

0.139 

1.18 


120 

0.193 

2.07 ± 0.50 

0.152 

1.27 


136 

0.212 

2.33 ± 0.45 

0.171 

1.24 

1300 

113 

0.211 

2.47 ± 0.27 

0.191 

1.11 


112 

0.214 

2.82 ± 0.41 

0.218 

0.98 


117 

0.189 

2.47 ± 0.46 

0.191 

0.99 

1400 

113 

0.201 

2.61 ± 0.29 

0.202 

1.00 


112 

0.191 

2.56 ± 0.19 

0.198 

0.96 


115 

0.199 

2.64 ± 0.57 

0.204 

0.98 

1200 

115 

0.128 

Dry O, 5 
1.76 ± 0.45 

0.129 

0.99 

1300 

112 

0.152 

1.83 ± 0.34 

0.141 

1.08 

1400 

112 

0.163 

1.87 ± 0.30 

0.145 

1.12 


♦Results are from work conducted in Ref. 10. Conversion from x to Awt assumes density of amorphous silica at 1200°C (2.2 
g/cm') and density of cristobalite at 1300° and 1400°C (2.32 g/cm'). 'Uncertainties reported for oxide thicknesses are la for at 
least 20 measurements at 1 mm spacing along SiC/oxide cross section. ‘Only data obtained in fused quartz reaction tubes 


are shown. 


Table III. Parabolic Oxidation Rate Constant, * p , Determined from the Paralinear 
Oxidation of SiC in 50% H 2 O/50% 0 2 Compared to Other Studies 


Ap (mg 2 /(cm 4 h)) 

Temp 

(°C) 

50% H 2 O/50% 0>, 
(this study) 

50% HjO/50% 0 2> 
(Ref. 12) 

5 X jt r (10% H-,0/90% 0 2 ) 
(Ref. i0) 


1200 

(3.80 ± 0.06) X 10 -* 
(2.72 ± 0.02) x 10 3 

1.14 X 10" 3 

1.27 X 10" 3 

1.28 X 10' 3 
1.05 X 10~ 3 

1.46 X 10 3 
1.55 X 10 3 
1.65 X 10‘ 3 


1300 

(1.37 ± 0.02) X 10 3 
(1.57 ± 0.03) X 10 3 
(4.21 ± 0.02) X 10 5 

1.85 X 10" 3 
1.72 X I0' 3 
2.08 X 10' 3 

1.98 X 10 3 
2.05 X 10 3 
1.53 x 10 3 


1400 

(3.93 ± 0.06) X 10” 3 
(5.57 ± 0.04) X 10' 3 

2.05 X 10” 3 
2.29 X 10” 3 
1.86 X 10" 3 

1.79 X 10 3 
1.63 X 10 3 
1.73 X 10 3 



A sample TGA curve for each type of experiment as well as 
the fitted curves are shown in Figs. 4 through 6. The linear rate 
constants determined by all three techniques are summarized in 
Table IV and plotted in Fig. 7. The linear volatilization con- 
stants measured for the fused quartz coupons and determined 
from the fit of the SiC paralinear data to the model equation are 
in terms of mg Si0 2 /(cm’-h), whereas the linear volatilization 


constant measured for the preoxidized SiC coupons is actually 
in mg SiC/(cm 2 -h), since the oxide thickness remains nearly 
constant. In Fig. 7, the linear volatilization constants from 
the fused quartz coupons and the paralinear fits have been 
multiplied by MW sic /MW si0! so that all volatility rates are 
expressed in mg SiC/(cm 2 -h) and a direct comparison between 
the three techniques can be made. While there is some scatter, 
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Fig. 3. Summary of parabolic oxide growth rates experimentally 
determined from the paralinear oxidation of SiC. The solid line is the 
linear regression result for the experimental results of this study only. 
Data from Refs. 10 and 12 are shown for comparison. 
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Fig. 5. Linear volatility of bulk silica at 1 200°C in 50% H 2 O/50% O z 
flowing at 4.4 cm/s. The solid line shows the volatility rate determined 
by linear least-squares analysis of the experimental data. 


the three techniques give similar results for the linear volatiliza- 
tion of SiC in water vapor. The mass flux of SiC has also been 
calculated as a function of temperature between 1200° and 
1400°C using Eq. (13), assuming Si(OH) 4 (g) is the volatile 
species and using vapor pressures calculated from the equilib- 
rium constant for Si(OH) 4 formation from Hashimoto 17 to deter- 
mine p. In both Eqs. (13) and (14), J = k x (MW Sl(0H)4 /MW Sj02 ). 
There is good agreement between the experimentally deter- 
mined values and those calculated based on Hashimoto ’s data. 
These types of calculations are often accurate only to within an 
order of magnitude. 

The temperature dependence of the linear volatilization con- 
stant is very low. An activation energy of 3 ± 75 kJ/mol (95% 
confidence interval) was obtained in this study. Though the 
uncertainty is large, this result is not statistically significantly 
different from the value of 57 ±2 kJ/mol reported by Hashi- 
moto 17 for the enthalpy of formation of Si(OH) 4 . 

The thickness of the oxide on preoxidized SiC coupons was 
measured after the second exposure in the 50% water vapor- 
containing environment. Typically, the oxide layer grew in thick- 
ness by an average of 9% (range = “18% to +48%) in this 
100-hour exposure, while a weight loss was always observed, 
averaging 0.25 mg/cm 2 (range = 0.13 to 0.36 mg/cm 2 ). The 
scatter in weight loss and scale thickness was independent 
of temperature. 



Fig. 4. Experimentally determined paralinear weight change kinetics 
for as-received SiC at 1200°C in 50% H 2 O/50% 0 2 flowing at 4,4 cm/s. 
The solid curve is the sum of Aw, and Aw 2 given by Eqs. (7) and (8) 
for k p and determined by nonlinear least-squares analysis of the 
experimental data. 



Time (h) 


Fig. 6. Linear volatility of preoxidized SiC at 1200°C in 50% H 2 0/ 
50% 0 2 flowing at 4.4 cm/s. The initial oxide thickness for this test 
coupon was about 8 pm. The solid line shows the volatility rate 
determined by linear least-squares analysis of the experimental data 
after 30 h. 
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Fig. 7. Summary of linear volatilization rates for SiC in 50% H 2 0/ 
50% 0 2 at 4.4 cm/s. The solid lines are calculated from Eqs. (13) 
and (14) for the experimental conditions. The dotted line is the linear 
regression result for the experimental data (all three methods). 
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Table IV. Experimentally Determined Values for the Linear Volatilization Rate Constant, 
*„ in 50% H 2 O/50% 0 2 at 4.4 cm/s 


Temp 

(°C) 

Paralinear oxidation of SiC 
(mg Si0 2 /(cm 2 -h)) 

Linear volatility of 
fused quartz 
(mg Si0 2 /(cm 2 *h)) 

Linear volatility of 
preoxidized SiC 
(mg SiC/(cm 2 -h)) 

1200 

(5.05 ± 0.06) X 10- 3 
(5.27 ± 0.03) X 10- 3 

1.69 x lO' 3 
1.55 x 10“ 3 

5.43 X 10 3 
2.16 x 10 3 

1300 

(1.60 ± 0.04) X lO' 3 
(2.68 ± 0.04) X lO' 3 
(5.07 ± 0.02) X lO' 3 

1.97 X 10 3 
2.07 x 10" 3 
1.47 X 10 3 

3.33 X 10~ 3 
4.13 x 10“ 3 

1400 

(6.41 ± 0.06) X lO' 3 
(7.75 ± 0.04) X 10~ 3 

2.00 X 10 3 
2.07 X 10~ 3 
1.43 X 10 3 

3.76 X lO’ 3 
3.32 X lO' 3 


(3) Scale Characterization 

The crystallinity of the oxide scale increased with tempera- 
ture as shown in Fig. 8. The scales at 1200°C were still com- 
pletely amorphous (within the sensitivity of XRD) even after 
100 h exposures in 50% water vapor. While this temperature 
dependence was expected, a surprising result was that the 
amount of crystallinity was no different than that found in dry 
oxygen under the same conditions. 30 There are many references 
in the literature that state that water vapor increases the crystal- 
lization rate of silica (summarized in Ref. 12). 

Typical SEM surface morphologies are shown in Fig. 9. 
Some bubbles are visible in the scale at 1200°C. Note the 
absence of cracks in the scale, indicating the scale is amorphous 
in agreement with the XRD results. The spherulitic structure 
and cracks in the scale progressively increased with tempera- 
ture at 1300° and 1400°C. Cracks in the scale are due to the 
(3-to-a transformation of cristobalite on cooling. While bubbles 
are no longer found, irregularly shaped pores are now observed. 

IV. Discussion 

(1) Scale Growth 

The value for the parabolic oxidation rate constant deter- 
mined from the paralinear equation represents the oxide growth 
rate that would be observed in the absence of silica volatility at 
that same temperature and water vapor partial pressure. The 
results pertaining to the parabolic rate constant as well as the 
scale morphology will be discussed in full elsewhere. 12 Prelimi- 
nary results indicate that the parabolic rate constant is approxi- 
mately proportional to P(H 2 0), indicating that scale growth in 
the absence of volatility is controlled by molecular diffusion of 
water vapor in silica. 11 The amount of crystalline phase appeared 
to be independent of the presence of water vapor. The origin of 
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Fig. 8. X -ray diffraction spectra for SiC oxidized for 100 h in 
50% H 2 O/50% 0 2 as a function of temperature. 


bubbles in the amorphous scale may be due to the increased 
pressures of gaseous reaction products formed by the accel- 
erated oxidation reaction in water vapor. These bubbles were 
observed only in amorphous scales, probably due to the decreased 
viscosity of amorphous silica containing high amounts of 
water. 31 These comments on the bubbles are only speculation 
at this point, and further discussion will be left to a future 
publication. 12 

(2) Paralinear Oxidation 

The first indication that paralinear kinetics are followed for 
SiC in 50% H 2 0 vapor came from the discrepancy between the 
measured Awt and the Awt predicted from the measured oxide 
thickness. In all the 50% water vapor exposures, there is not 
enough weight gain to account for the amount of oxide found 
on the surface, even though the loss of C as CO is already 
accounted for. Therefore, some weight loss must be due to loss 
of silicon. Additional evidence for the volatility of the silica 
scale is the deposition of silicon- and oxygen-containing mate- 
rial on cooler surfaces in the reaction tube and on the sample 
hanger. The shape of the 100 h TGA curves themselves, while 
not exhibiting the long-term linear weight losses characteristic 
of paralinear behavior, do show the initial stages of paralinear 
kinetics: a significant leveling off of the weight gain or a maxi- 
mum in the weight gain is observed after 30 to 40 h. Experi- 
ments to times approaching the limiting oxide thickness (on the 
order of 1000 h) were not possible, due to limitations in the 
water vapor saturation system. Proof of silica scale volatility 
comes from the linear weight losses observed for preoxidized 
SiC with thick silica scales and the agreement of these rates 
with those found for volatility of bulk silica. 

The question may be asked why wasn’t any evidence of 
this paralinear behavior previously observed for silica-forming 
materials? Much of the previous work (summarized in Ref. 10) 
was conducted for short times, less than 30 h, in which the 
linear volatilization is overwhelmed by the parabolic oxidation 
of SiC. Secondly, in previous studies, impurity effects were 
not isolated from water vapor effects. Impurities enhance the 
parabolic oxidation rate of SiC, 30 but are expected to have 
little effect on a boundary-controlled volatilization rate. By 
examining Eqs. (7) and (8), it can be seen that an increase in k p 
with a constant £, would increase the time at which the maxi- 
mum in weight gain occurs. Volatility effects on the weight 
change would be more difficult to observe under these condi- 
tions. Finally, many of the earlier investigations were also 
conducted in low partial pressures of water vapor («*3% to 
10% H 2 0) at low gas velocities (^0.4 cm/s) where volatility is 
not observed. 

However, there were some hints that volatility of the silica 
scale would be important under the conditions studied here. 
Ligenza 20 observed the solubility of the oxide formed on silicon 
in the stream under a variety of conditions, for example, 650°C, 
75 atm, 2 h. Deal 21 found that in wet oxygen (85% H 2 0) and 
steam at 1 atm total pressure the amount of silica formed on 
silicon after 4 h was less than that expected from the amount of 
silicon consumed. This was attributed to the erosion of the 
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Fig. 9. Scanning electron micrographs showing the oxide surface morphology of CVD SiC oxidized in 50% H 2 O/50% 0 2 for 100 h as function of 
temperature: (a) 1200°C, (b) 1300°C, and (c) 1400°C. 


oxide in the wet oxidizing atmospheres. The mechanism of 
erosion was not determined, nor was the rate of erosion quanti- 
fied. Choi et al. n observed that “parabolic rate constants” for 
Si 3 N 4 calculated from oxide thicknesses after exposure in water 
vapor at 1 100° to 1200°C for times around 10 h and less actually 
decreased with water contents over 40% to 60%. This behavior 
may be due to paralinear oxidation rather than the blocking 
effect of the NO or NH 3 reaction products proposed by the 
authors. Silicon nitride may show paralinear behavior at shorter 
times, since the parabolic oxidation rate is much lower than that 
of SiC. Paralinear kinetics for Si 3 N 4 have not yet been studied. 

Since the conditions of this study most closely match the 
conditions studied by Hashimoto 17 (i.e., 1 atm total pressure 
H 2 0/0 2 mixtures, P(U 2 0) varying between 0.1 and 0.6 atm, 
T = 1100° to 1500°C), it was assumed throughout this study 
that Si(OH) 4 is the major volatile species. Preliminary mass 
spectrometric results obtained for Si0 2 at 1 atm total pressure, 
90% to 95% H 2 0 in 0 2 , and 1300°C confirm the presence of 
this vapor species. 18 As mentioned in the Results section, there 
is relatively good agreement between the linear volatilization 
constants determined by the three experimental methods and 
the values calculated for Si(OH) 4 volatility using Hashimoto’s 
data. This agreement supports the assumption that Si(OH) 4 is 
the volatile species formed under the conditions of this study. 

The formation of volatile species from the fused quartz reac- 
tion tube apparently did not saturate the gas stream. If the fused 
quartz tubes at least partially saturated the gas stream and 
suppressed volatilization of the test coupon, then the experi- 
mental linear weight loss rates would be less than those calcu- 
lated based on the boundary layer equation under nonsaturated 
conditions. 

The nonlinearity in weight loss observed at the beginning of 
the exposures of the preoxidized SiC coupons may be due to 
sealing of cracks and imperfections in the thick oxide scale 
which were observed by SEM after preoxidation. After irregu- 
larities in the scale are healed, only linear volatility would be 
observed. A limiting oxide thickness calculated from all k p and 
k { results at 1200°C in Tables III and IV give a limiting oxide 
thickness of 6.9 |xm. Preoxidized samples had oxidized thick- 
nesses of 6 to 8 pm, which is close to the limiting oxide 
thickness expected. Therefore, 9% growth in oxide thicknesses 
during subsequent 100 h exposures in H 2 0/0 2 indicates growth 
is minimal, as would be expected for scales near the limiting 
oxide thickness. 

The volatilization rates and activation energies did not 
depend on the crystallinity of the silica scale. The volatility 


rates of two fused quartz coupons were measured, one in an 
alumina tube (coupon crystallizes in a few hours) and in a fused 
quartz tube (coupon remains amorphous) at 1200°C for 100 h. 
The volatility rates were within a factor of 2 of each other. This 
is the same magnitude of variability found in exact repeat experi- 
ments. Even if the reaction rates for water with amorphous and 
crystalline silica differed, no difference in volatility rates would 
be expected for a boundary layer-controlled process. 

(3) Combustor Applications and Life Predictions 

The implications of this study for the long-term chemical 
durability of SiC in combustion environments will now be 
addressed. First, of most importance, is the conclusion that the 
volatility of silica controls the lifetime of any SiC component. 
At long times (on the order of hundred hours), a steady state 
should be achieved in which a constant thickness of oxide exists 
on the SiC. The oxidation of SiC then proceeds linearly at the 
same rate that the silica volatilizes. The SiC surface recession 
also occurs at a linear rate. The only influence of the parabolic 
oxidation rate constant, k p , on this behavior is to establish the 
limiting thickness of the oxide layer which recedes into the 
SiC. For long-term applications, the recession rate of SiC can 
therefore be approximated using the linear rate constant alone. 
This has also been demonstrated for other paralinear processes 
such as alumina formation and spallation on NiAl. 33 

The effect of the relative temperature dependence of the rate 
constants on scale growth should also be considered. Because 
the activation energy for oxidation is greater than that for vola- 
tilization, the limiting oxide thickness should increase with 
temperature. An additional point is that silica volatility will be 
a problem even at the lower application temperatures of SiC 
(1200°C and below). As the temperature is lowered, the oxida- 
tion rate would decrease compared to the volatility rate. At 
some point, paralinear kinetics would no longer apply. SiC 
oxidation would be dominated by linear oxidation kinetics 
due to the thinness of the oxide layers formed. The oxidation/ 
volatilization reaction kinetics under these conditions would be 
described as the summation of two linear processes. 6 

Given a recession rate goal, extrapolations of the results from 
this paper can be made to predict lifetimes of SiC in real 
applications. The paralinear behavior observed at short times is 
assumed to be negligible for the long-term lifetime of SiC in 
water vapor-containing environments. Thus, Eq. (13) can be used 
to predict the SiC recession rate for long times. Equation (13) 
can also be reduced to a simple relationship between the mass 
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flux, k u and process variables such as water vapor partial pres- 
sure, total pressure, and gas velocity. The temperature depen- 
dence is not included in this relationship, first, because the 
results from this study were conducted at typical desired use 
temperatures of SiC so no temperature extrapolations are 
needed, and second, the temperature dependence was found to 
be minimal. The velocity dependence of the flux is clearly v 1/2 ; 
however, the pressure dependence is more complex. The diffu- 
sion coefficient, D, is proportional to 1 /P toMt p' a P tmii , and, 
assuming Si(OH) 4 is the volatile species, according to Eq. (4), 
p P(H 2 0) 2 . The resulting simplified relationship is 


v 1/2 P( H 2 0) 2 

^1 a p 1/2 

r total 


(17) 


Under combustion conditions, for example: T ~ 1200°C, 
P tolal = 10 atm, P(H 2 0) = 1 atm, and gas velocity = 3000 cm/s, 
the volatility rate would be about 33 times greater than deter- 
mined in this study. In addition, these calculations were based 
on laminar flow conditions limited by diffusion through a gas 
boundary layer. The Reynolds number (for flow in a smooth 
pipe geometry) for the combustion conditions above is about 
3 X 10\ which is greater than 2100 where the transition to 
turbulent flow begins. 34 Therefore, the actual combustion condi- 
tions may be turbulent. The estimated volatility rates based on 
Eqs. (13) or (17) serve as a lower bound. The volatility rate 
calculated from the Langmuir equation, Eq. (14), for vapori- 
zation in a perfect vacuum can serve as an upper bound. 
For example, this upper bound is shown in Fig. 7 for the 
1 atm conditions studied here. Actual volatilization rates will 
fall between these two bounds, depending on the turbulent 
conditions. 

Mass fluxes can easily be converted to SiC recession rates 
using the density of SiC (^3.2 g/cm 3 ). Using the average k i 
from this study from all temperatures between 1200° and 
1400°C, about 2.5 X I0~ 3 mg/(cm 2 -h), extrapolated to the com- 
bustion conditions above, a recession rate of 2.6 X 10 5 cm/h is 
calculated, assuming laminar flow. Recession rates would be 
greater in the case of turbulent flow. Lifetimes of SiC compo- 
nents in combustion environments can then be computed based 
on the maximum acceptable amount of component recession. 
The volatility of silica scales on SiC in combustion environ- 
ments is clearly a concern for thin components. 


V. Conclusions 

It has been demonstrated that SiC undergoes paralinear 
oxidation/volatilization reactions in an environment containing 
50% H 2 O/50% oxygen at 1 atm total pressure and temperatures 
between 1200° and 1400°C. Based on this work, the results of 
Hashimoto, 17 and preliminary atmospheric mass spectrometry 
results, 18 it is believed that the silica scale volatilizes as an 
Si(OH) 4 (g) species. Because the volatility of silica in water 
vapor controls the recession rate of SiC, lifetimes of SiC based 
on a component thickness criterion can be predicted based on 
the linear volatilization rate of Si0 2 alone. Extrapolations to 
hypothetical combustion conditions have shown that the life- 
time of thin SiC components in some combustion environments 
are inadequate for long-term applications. 
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